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Abstract. We use a smoothed version of the explicit formula to find an approximation to the 
Riemann zeta function as a product over its nontrivial zeros multiplied by a product over the 
primes. We model the first product by characteristic polynomials of random matrices. This provides 
a statistical model of the zeta function that involves the primes in a natural way. We then employ 
the model in a heuristic calculation of the moments of the modulus of the zeta function on the critical 
line. This calculation illuminates recent conjectures for these moments based on connections with 
random matrix theory. 



1. Introduction 

An important theme in the study of the Riemann zeta function, C( s )> nas been the estimation 
of the mean values (or moments) 



4(T) = iy |C(^ + ii)| 2 " dt. 



These have applications to bounding the order of £(s) in the critical strip as well as to estimating 
the possible number of zeros of the zeta function off the critical line. Moreover, the techniques 
developed in these problems, in addition to being interesting in their own right, have been used 
to estimate mean values of other important functions in analytic number theory, such as Dirichlet 
polynomials. 

In 1918 Hardy and Littlewood |Hj proved that 

/iCn-iogT 

as T — > oo. Eight years later, in 1926, Ingham JO] showed that 

I 2 (T) ~ JL (log Tf. 

There are no proven asymptotic results for when k > 2, although it has long been conjectured 
that 

I k {T) ~ c k {\ogTf 

for some positive constant Cfe. Conrey and Ghosh (unpublished) cast this in a more precise form, 
namely, 



where 



«T)~«^ogr)<' 



1 \ v-^ / 1 [m + k) 
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the product being taken over all prime numbers, and g(k) is an integer when k is an integer. The 
results of Hardy-Littlewood and Ingham give g(l) = 1 and g(2) = 2, respectively. However, until 
recently, no one had formed a plausible conjecture for g{k) when k > 2. Then, in the early 1990's, 
Conrey and Ghosh [Ij conjectured that g(3) = 42. Later, Conrey and Gonek jS] conjectured that 
g(4) = 24024. The method employed by the last two authors reproduced the previous values of 
g{k) as well, but it did not produce a value for g(k) when k > 4. 

It was recently suggested by Keating and Snaith ^3] that the characteristic polynomial of a large 
random unitary matrix can be used to model the value distribution of the Riemann zeta function 
near a large height T. Their idea was that because the zeta function is analytic away from the 
point s = 1, it can be approximated at s = ^ + it by polynomials whose zeros are the same as 
the zeros of £(s) close to t. These zeros (suitably renormalized) are believed to be distributed 
like the eigenangles of unitary matrices chosen with Haar measure, so they used the characteristic 
polynomial 

N 

Z N (U,9) = l[(l-e^^), (2) 

n=l 

where the 9 n are the eigenangles of a random N x N unitary matrix U, to model C( s )- For scaling 
reasons they used matrices of size N = logT to model £(^ + ti) when t is near T. They then 
calculated the moments of \Zn(U, 9)\ and found that 

e n [\z N (u,e)\ 2k ] ~ H^y^ 2 > (3) 

where E^r denotes expectation with respect to Haar measure, and G(z) is Barnes' G-function. 
When k = 1,2,3,4 they observed that 

G 2 (k + 1) _ g(k) 
G{2k + 1) r(A: 2 + l)' 

where g(k) is the same as in the results of Hardy-Littlewood and Ingham, and in the conjectures of 
Conrey-Ghosh and Conrey-Gonek given above. They then conjectured that this holds in general. 
That is, they asserted 

Conjecture 1 (Keating and Snaith). For k fixed with Refc > —1/2, 

as T — » oo, where a(k) is given by (^Q) and G is the Barnes G-function. 

The characteristic polynomial approach has been successful in providing insight into other im- 
portant and previously intractable problems in number theory as well (see, for example, j!6j for a 
survey of recent results). However, the model has the drawback that it contains no arithmetical 
information — the prime numbers never appear. Indeed, they must be inserted in an ad hoc man- 
ner. This is reflected, for example, by the absence of the arithmetical factor a(k) in equation ©. 
Fortunately, in the moment problem it was only the factor g(A;), and not a(k), that proved elusive. 
A realistic model for the zeta function (and other L-functions) clearly should include the primes. 

In this paper we present a new model for the zeta function that overcomes this difficulty in 
a natural way. Our starting point is an explicit formula connecting the zeros and the primes 
from which we deduce a representation of the zeta function as a partial Euler product times a 
partial Hadamard product. Making certain assumptions about how these products behave, we 
then reproduce Conjecture ^ Our model is based on the following representation of the zeta 
function. 
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Theorem 1. Let s = a + it with a ^ and \t\ ^2, let X ^ 2 be a real parameter, and let K be any 
fixed positive integer. Let u{x) be a nonnegative C°° function of mass 1, supported on [e l ~ l l x ,e], 
and set 



OC 



U(z) = I u(x)E 1 (zlogx) dx , (4) 
Jo 

where E\(z) is the exponential integral e~ w /w dw. Then 

C(s) = P x (s)Z x (s) [ 1 + ^ l ^ )K ^+0(X-^logX^j , (5) 



where 



n s log n 



A(n) is von Mangoldt's function, and 

Z x (s) = exp(-Y / U[(s-p n )logX)) . (7) 

V Pn I 

The constants implied by the O terms depend only on u and K. 



We remark that Theorem ^ is unconditional — it does not depend on the assumption of any 
unproved hypothesis. Moreover, it can easily be modified to accommodate weight functions u 
supported on the larger interval [1, e\. Finally, as will be apparent from the proof, the second error 
term can be deleted if we replace Px(s) by 

P X ( S ) = exp ( y J^L„ (e iog«/iogX; 
w \ ^ n s \ogn K ' 

where v(t) = u(x)dx. 

To clarify we temporarily assume the Riemann Hypothesis (RH) and take s = h + it. We 
shall denote the nontrivial zeros of ((s) by p n = | + i7n, ordered by their height above the real 
axis, with j_ n = —j n . Since the support of u is concentrated near e, U(z) is roughly E\(z), which 
is asymptotic to —7 — logz as z — > 0. Here 7 = 0.5772... is Euler's constant. Thus, for those 
ordinates 7 n close enough to t, we see that 

exp ( - U(i(t - 7n ) log X)) « 1 (t - 7n ) log X . 

We expect the ordinates farther away not to contribute substantially to the exponential defining 
Zx(s). Now, Px( s ) ra I1p<a(^ — P hence our formula looks roughly like 

C(i+>*)« n^-P"*" 1 *) -1 II (l(*-7n)^lo g x]. (8) 

[t- 7n [<l/logX 

This formula is a "hybrid" consisting of a truncated Euler product and (essentially) a truncated 
Hadamard product, with the parameter X mediating between them. Near height T we are approx- 
imating part of the zeta function by a polynomial of degree about log Tf log X. The rest of the 
zeta function, which comes from the zeros we have neglected, is approximated by the finite Euler 
product. Formally, when we take X large, we reduce the number of zeros used to approximate 
zeta, but make up for it with more primes; and when we take X small, we approach the previous 
model (|2"|). Note however, that in order for the error terms in (JSJ) to be smaller than the main term, 
it is necessary to work in an intermediate regime, where both the zeros and the primes contribute. 
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To see how to use the model, and as a test case, we heuristically calculate Ik(T). The new 
model is more elaborate than the original one, so more work is required. Nevertheless, the idea 
is straightforward. The 2/cth moment of \C(h + ^)l is asymptotic to the 2/cth moment of l-Px^ + 
it) Zx{\ + rt) | - We argue that when X is not too large relative to T, the 2/cth moment of this 
product splits as the product of the moments. We call this the "Splitting Conjecture". 

Conjecture 2. (Splitting Conjecture.) Let X and T — > oo with X = 0((logT) 2 ~ e ). Then for 
k > —1/2 we have 

1 /-2T / i r 2T \ / i r 2T 



y , |C(s + it)| dt ^[fJ T \ Px ^ + ^\ dt ) x {fJ T \ Zx ^ + lt )\ dt 

In Section El we calculate the moments of P rigorously and establish the following theorem. 

Theorem 2. Let 1/2 ^ c < 1, e > 0, and let k be any real number. Suppose that X and T — > oo 
andX = ((logT) 1 ^ 1 "^ 6 )) . T/ien we have 

1 f 2T / / 1 



y |Px(^ + it)| 2fc dt = a(A;, CT )F x (A:, C T) (^1 + 0* 



uniformly for c ^ a ^ 1, where 




m\2 



and 



p / „ P 

m=0 



Fx(k,cr) 



((2a) k e- k E ^ 2 °-V tos JO j/ CT > 1/2 , 
( e TlogX) fe2 if a = 1/2. 

Here E\ is the exponential integral, and 7 = 0.5772 . . . is Euler's constant. 



Note that a(fc, 5) is the same as a(k) in (|T]). 

In Section0]we conjecture an asymptotic estimate for J^ T \Zx{\ + it)\ 2k dt using random matrix 
theory. We introduce random matrix theory in the following way. The statistical distribution of 
the ordinates 7™ is conjectured to coincide with that of the eigenangles 8 n of N x N random unitary 
matrices chosen with Haar measure for some N (see for example |17j . |19j and |15|). The choice 
of N requires consideration. The numbers 7„ are spaced 27r/logT apart on average, whereas the 
average spacing of the 9 n is 2n/N, and so we take N to be the greatest integer less than or equal to 
logT. We therefore conjecture that the 2kth moment of \Zx{\ + when averaged over t around 
T, is asymptotically the same as \Zx{\ + 1 t)\ 2k when the 7„ are replaced by 9 n and averaged over 
all unitary matrices with N as specified above. We perform this random matrix calculation in 
section |1J and so obtain the following conjecture: 

Conjecture 3. Suppose X, T — » 00 with X = 0((logT) 2 ~ e ). Then for any fixed k > —1/2, we 
have 

dt G 2 (k + 1) ( logT ^ 2 



T J T 1 Vi " G(2k + 1) Ve^logX 

We actually expect conjecture 131 to hold for a much larger range of X, but the correct bound on 
the size of X with respect to T is unclear. 

We note that this asymptotic formula coincides with that in (j3J) when there A?" is taken to be 
on the order of logT/e 7 logX. This is consistent with the fact that the polynomial in (JHJ) is of 
about this degree. Alternatively, the mean density of eigenvalues is N divided by 2ir, and this is 
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comparable to the mean density of the ordinates of the zeros when multiplied by e^logX, as they 
are in (jHJ). 

Combining the result of Theorem |^1 with the formula in Conjecture |3] and using the Splitting 
Conjecture, we recover precisely the conjecture put forward by Keating and Snaith. Note that, as 
must be the case, all X-dependent terms cancel out. 

In Section [SJ we prove 

Theorem 3. Let e > and let X and T -> oo with X = 0((logT) 2_e ). Then for k = 1 and k = 2 



Since £(± + d)Px(\ + i*) -1 = Z x {\ + it) (1 + o(l)) for t e [T,2T], it follows from this that 
Conjecture |31 holds when k = 1 and k = 2. Moreover, combining Theorem |3] with our estimate for 



from Theorem |21 we also see that Conjecture El holds for k = 1 and k = 2. Thus, we obtain the 
Corollary. Conjectures^ and\^ are true for k = 1 and k = 2. 

Clearly our model can be adapted straightforwardly to other L-functions (see |14| ). It can also 
be used to reproduce other moment results and conjectures, such as those given by Gonek jS] and 
by Hughes, Keating and O'Connell jHl concerning derivatives of the Riemann zeta function at the 
zeros of the zeta function. We also expect it to provide further insight into the connection between 
prime numbers and the zeros of the zeta function. It would be particularly interesting to determine 
whether the model can be extended to capture lower order terms in the asymptotic expansions of 
the moments of C(l/2 + it) and other L-functions, c.f. 0. 



We begin the proof by stating a smoothed form of the explicit formula due to Bombieri and 



Lemma 1. Let u(x) be a real, nonnegative, C°° function with compact support in [l,e], and let u 
be normalized so that if 



we have 





2. The Proof of Theorem 1 



Hejhal 0. 




then v(0) = 1. Let 



Jo 

be the Mellin transform of u. Then for s not a zero or pole of the zeta function, we have 





m=l 



(10) 



where the sum over p runs over all the nontrivial zeros of the zeta function. 
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This lemma is proved in a familiar way, beginning with the integral 

i r C dz 

— U z + s )u{l + z\ogX)-, 
2vriy (c) C z 

where the integral is over the vertical line Hez = c = max{2, 2 — Re s}. 

The support condition on u implies that v(e logn ^° sX ) = when n > X, so the sum over n is 
finite. Furthermore, if |Imz| > 2, say, then integrating u by parts K times, we see that 



\u{z)\ ^max|«^(x)| 



T(z) 



T(z + K) 



, e Rcz+K + ^ 



e max{Rcz+K,0} 

^ max |it^^(x)| 



(i + N) 



K 



for any positive integer K. Thus, the sums over p and m on the right-hand side of (jlUj) converge 
absolutely so long as s ^ p and s ^ —2m. This, in fact, is the reason we require smoothing. 

Next we integrate (fTU)l along the horizontal line from so = <7o + to +oo, where cjq ^ and 
|io| ^ 2. If the line does not pass through a zero, then on the left-hand side we obtain — logC(so). 
We choose the branch of the logarithm here so that linio-^oo log£(s) = 0. If the line of integration 
does pass through a zero, we define log £(<7 + ii) = lim e ^ + ^ (log£(cr + i(t + e) + logger + i(t — e)) . 
Recalling the definition of U(z) in Q), we see that 

f°° u(l - (s-z) log X) r ( \w(( \ l Yl ^ 

/ as = / ulx)h/i{{sQ — z) log X logx) ax 

J s s ~ z Jo (12) 

= U((s -z) log X), 

provided that so — z is not real and negative (so as to avoid the branch cut of E\). If it is, we use 
the convention that U({s$ — z) log X) = lim e _ >0 + \ (U((sq — z) log X + le) + U ((sq — z) log X — ie)) . 
Note that the logarithms in (|12j) are both positive since the support of u is in [l,e] and X ^ 2. It 
therefore follows from that 



logC(^o) = £ J^-«(e logr,/logl ) - £ I7(( i0 - p) logX) 
z — ' n s ° log n 



(13) 

CO v ' 

+ U((s - 1) log X)-J2 U((s + 2m) log X) . 

m=l 

The interchange of summation and integration in the sums is justified by absolute convergence. 
This representation holds for all points in Re s ^ not equal to the pole or one of the zeros of the 
zeta function. 

We next suppose that the support of u is contained in [e 1_1// ^, e] with the same X as in (|13|) . It 
is easy to see that there is a smooth nonnegative function / with support in [0, 1] and total mass 
one such that u(x) = Xf(X\og(x/e) + l)/x. Since max^ \ f( K \x)\ is bounded and independent of 
X, we see that max^ |u^(x)| <C_k" X k+1 . It therefore follows from that 

gmax{o-,0} 

=00 «* (1 + M) * ■ 
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Prom this and (|12|). and since |£q| ^ 2, we find that if r is real, then 

tt« \ i in /•°° n(l-( g -r)logX) ^ 

J7((s - r) log X) = / ds 

J s s r 

f 00 ^max{r-(j, 0} 

<<K JteXY? J ao \(a - r) + it \K+i da 

j^-_ft"+l+max{r— (jo, 0} rco ^ 
<K 7T— / 77- ,v , r: fyZi dcJ 



(log*)* 7 CT0 Ka-O+rfol^ 1 

j^-fC+l+max-jY— Co, 0} 



(|a -r|logX)* ' 
In particular, for any fixed positive integer ii" we have that 

j^_ftT+l+max{l— Co, 0} 

C(0»o-l)log*)<* (|8o|log ^ 

and, since do ^ 0, that 

X K+1 



X + 1 

S U « S0 + 2m) bg X) <<K (bgTF ^ | 80 + 2m|^i 

m=l m=l 1 



(| So |logX)* • 

Inserting these estimates into 1)13(1 and replacing so by s, we find that 
for cr^O, |i| ^ 2, and if any fixed positive integer. Exponentiating both sides, we obtain 

x K+2 



C(s) = P x (s)Z x {s) (l + O 

where 



\s\\ogX) K 



P x (s) expf £ -A^L- i; . 



n s log n 



and 



Z x (s) = exp ^- ^ U((s - p) log X) j 



We now wish to show that replacing Px(s) by 



Px(s) = exp 



A(n) 



ra s log n 
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only introduces a small error term into (|14l) . To see this, note that w (( e iogn/lo g x^ = i f or n ^ 
X^ l / X because the support of u(x) is in [e l ~ l l x , e\. Therefore, 

A(re) / i ogX/logn , _ 
t—' n s log n V 




E : 

<exp(X- CT logX) . 

This completes the proof of Theorem ^ provided that s is not a nontrivial zero of the zeta 
function. To remove this restriction, we recall the formula 



m! m 

m=l 



where | argz| < 7r, logz denotes the principal branch of the logarithm, and 7 is Euler's constant. 
From this and @ we observe that we may interpret exp(— U(z)) to be asymptotic to Cz for some 
constant C as z — > 0. Thus, both sides of © vanish at the zeros. 



3. The Proof of Theorem [21 



We begin with several lemmas. 
Lemma 2. Let X > 2 and se£ 



*M-nH)"' n 



Then for k any real number we have 

P x (s) k = P x {s) k (l + O k f 1 

uniformly for a ^ 1/2. 
Proof. By @ we have 



logX 



\ & / ps?A \ l^j^N p J ^ / 

where A/ p = [log X/ log p], the integer part of \ogX/\ogp. Therefore 

The primes < p ^ X have iV p = 1, and we note that the j = 2 term for these primes in the 
first double sum exactly cancels the j = 1 term in the second. Hence the argument in the exponent 
is 

^ W I p a(N p + l) + 
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Now p N p +1 > X since N p + 1 > log X/ log p, so, for a ^ 1/2, this is 



+ 



It follows that 



as required. 



logX X 1 / 4 logX J logX' 
Px (s) k P* x (s)- k = l + O k , losX 



Lemma 3. Let k be a real number. Let 1/2 ^ c < 1 be arbitrary but fixed, and suppose that 
2 ^ X <C (logT) 1 /( 1 ~ c+e ), where e > is also fixed. Then 



1 

T 



2T 



\P x (a + it)\ 2k dt = a(k,a) ]J ( 1 - -L 



-fe 2 



l + O fc 



logX 



uniformly for c ^ a ^ 1, where a(k, a) is given by ©• 



Proof. We write 



E 

n=l 



11(1 "f 



n a+^r 

\/A<p^X 



(15) 



Let <S(X) denote the set of X-smooth numbers, that is, S(X) = {n : p \ n ==>- p X}. Then 
a k{ n ) = d>k( n )i the feth divisor function, if n € <S(\/X); aifc(p) = d k {p) for all p ^ X; and afc(n) = 
if n <5(X). It is also easy to see that 



(i-p-r fc (i+K 2 r fc 

Comparing this with 



1 



1 1 + 



1A 1 



exp /c 1 =- + 

\ \p s 6p 6s 4p 4s 



fi.s 



we find that for k ^ 0, yX<p^X, and j = 1, 2, . . ., 

«S ot k (jp?) < d 3k/2 (p i ) , 

while for < 

|«fc(y)| < a|fe|(p J ") < d3|fc|/2(p f ) • 

We now truncate the sum in (|15|) at T e , where 9 is a small positive number to be chosen later, 
and obtain 



E 

ndS(X) 
n<T s 



a k (n) 



( 



+ 



E 

neS(X) 
- n>T e 



|a fc (n)| 



\ 



rr 
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For e > fixed and a ^ c, the sum in the O-term is 



« £ 



f n y d m/2 {n) <T _ e9 ^ d 3 \ k \ /2 {n) 



n>T B 



n" 



neS(X) 



T ~ ee IK 1 - p e ~T = T ' ee ex p \\ k \ E p e 

<T" ey exp(0' 



.(l-c + e)logX, 

Now suppose that 2 ^ X <C (logT) 1 ^ 1- ^) with the same e. Then this is 

■\k\ log T" N 



< T~ ew exp O 



log log r 



Thus, we find that 



P x (s) k 



Ea k (n) 
— — "M'A- 



+ o, (r-^/ 2 ) . 



ne5(X) 
n<T e 



(16) 



Next we calculate ^ J -2 |P£(s)| 2fc di. By Mont gomery and Vaughan's mean value theorem for 
Dirichlet polynomials .IB], we have 



2T 



neS(X) 



dt=(T + 0(T e )) 



ak{n) 

2a 



n4T e 
neS(X) 



n 



Using the method above, we may extend the sum on the right to infinity with an error again no 
larger than O k (T~ ee / 2 ). Thus, taking 6 = 1/2, say, we find that 



2T 



a k (n) 



neS(X) 

i r 2T 



n£5(X) 



(17) 



We next note that if = A J" T |aj(t)| 2 di , i = 1, 2, and A\ ^ 0, then 

^ £ T \ ai (t) + a 2 (t)\ 2 dt = A 1 (l + ((AsMO 1 / 2 ))) + A 2 , 

the O-term arising from the Cauchy-Schwarz inequality applied to the "cross term". We use this 
with a\(t) the sum on the right-hand side of (|16f) and a 2 (t) the error term (with 9 = 1/2). Since 
ccfc(l) = 1, we see from (I17j) that A% 3> 1. It therefore follows that 



i r2T 

f J T \P* x (a + it)\ 2k dt=(l + O k (T-^)) 



otk(nf 



neS(X) 



n 



2a 



(18) 
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Since ak(n) = dk(n) for n € S(y^X), and ajt(p) = dk{p) for y/X < p ^ X, we may write the sum 
as 

a t (n) 2 rr a t (j>">) 2 \ 



E =r-n E 



n£S(X) pS^A \m=0 F / 

rr /f^ 4(p m ) 2 \ n / 4(p) 2 f Wy m ) 2 \ 

11 \ 2^1 v 2ma I 11 + + Zj v 2ma I 

_ tt (f 4W 2 ) tt ((, , 4(p) 2 f Wp m ) 2 \ /f 4W 2 l 

1 1 I Z-/ B 2niff I 11 11 r) 2<T B 2m<T J / B 2m<T f ' 

p<;a \m=0 ^ / y/X<p^X \ V m = 2 / m =o J 

Factoring 1 + d) c (p) 2 /p 2eT (which is at least 1) out of the numerator and denominator of the last 
product, we see that the product equals 

n (i + o t (£)) = ex P ( 0k (*£.)) = i + o, 

Hence, 



— P / / \ \ log -X" / / V log X 

v / X<p<A \ \ / / V 



II (S ^Lr ) (i + Ofcf ! X 

neS(A) p^A \m=0 F / V V B / 

Writing the product here as 

-k 2 



HU 1 p^) E^jnf 1 p^) 



p<A \ v ^ 7 m=0 ^ / p^A 

we note that the first of the two factors may be extended over all the primes, because 

p>X \ V F 7 m=0 F / p>A V \f / / 



=1 + o k 



logX 

Thus, by the definition of a(k,a) in @, we find that 



logX 



ne-S(A) p^A 

The lemma follows from this and Q18|). 
Lemma 4. //A; is a real number, then 

n(i-^)""^x(^)(i +0t(i ^)) 

uniformly for a ^ 1/2, where 

„ n s r c(2(J )^ e -^i((2-i)iogA) if<T>1/2i 

[(e 7 logXf if a = 1/2, 

and Ei is the exponential integral. 
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Proof. Mertens' theorem asserts that 

-l 



nH)~ =^(1+0^)) 



p^X 

Raising both sides to the k 2 power establishes the result when a = 1/2. When a > 1/2, we see 
that 

By the prime number theorem in the form ip(x) = X^n<ir^( n ) = ^ + 0(x/(\ogx) ), we find that 



?M 1_ ^) = ~It(^ +0 (^)) 



p>X v ^ ' p>X 



00/1 ^ f 1 A A du 1 



x \u 2a \u 4a J J logu V( lo g^) A 
-Si((2a-l)logX) + 0' ' 



(logX)^ 

" =C(2afexp(-^ 1 ((2a-l)logX))(l + O fc ( i -^)) 



Hence, 

-it 2 



p^X 

as asserted. 

The proof of Theorem |21 now follows immediately from Lemmas El El and El 

4. Support for Conjecture El 
In this section we give heuristic arguments supporting Conjecture El which we restate as 

as T — ► oo, where Zx(s) is given by 0. 

We assume the Riemann Hypothesis. Since Re-E'i(ix) = — Ci(|x|) for where 

c&(*)=- r— ^, 

Jz W 

we find that 

1 C 2T 1 f' 2T ( f' e \ 

- J t \Z x {\ + it)\ 2k dt=- ]Je W (2k J u{y) Ci(|i- 7n | logy log X) dy\ dt , (19) 

where is a smooth nonnegative function, supported on [e 1_1// ^,e] and of total mass 1. Since 
the terms in the exponent decay as \j n — t\ increases, this product is effectively a local statistic. 
That is, the integrand depends only on those zeros close to t. In recent years considerable evidence 
has been amassed suggesting that the zeros of the Riemann zeta function around height T are 
distributed like the eigenangles of unitary matrices of size log T chosen with Haar measure (see, for 
example, the survey article ^SJ). We therefore model the right-hand side of (|T§|) by replacing the 
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ordinates 7 n by the eigenangles of an N x N unitary matrix and averaging over all such matrices 
with Haar measure, where ./V = [logT]. Thus, the right-hand side of l)19j) should be asymptotic to 

Y[exp(2k [ u(y) Ci(\9 n \ logy log X) dy 



E 



iV 



n=l 



where the 9 n are the eigenangles of the random matrix and E^r [•] denotes the expectation with 
respect to Haar measure. However, since the eigenangles of a unitary matrix are naturally 2ir- 
periodic objects, it is convenient to periodicize our function, which we do by defining 



(f)(9) = exp 2k / u(y) 



oo 

E 

j=-O0 



Ci(|0 + 27rj| logy log X) dy 



(20) 



It will follow from our proof of LemmaElthat the terms with j ^ 0, which make the random matrix 
calculation much easier, only contribute <Cfc 1/logX to (f)(9) when — tt < 9 ^ tt. Hence they do not 
affect the accuracy of the model. Thus, we argue that 

N 



- J t \z x (\+t) 



2k 



dt ~ E 



N 



n 

,n=l 



(21) 



The remainder of this section is devoted to the proof of 

Theorem 4. Let (f)(9) be defined as in (|2Uj). then for fixed k > — | and X ^ 2, we have as N — > oo, 

\2 / »r \ fc 2 



E 



iV 



■ N 

n 

.n=l 



N 



(G(k + 1)Y 
G(2k + 1) \e~<logX 



l + O h 



1 



log* 



Remark. The random matrix model of Keating and Snaith |13j for the moments of the Riemann 
zeta function involved the characteristic polynomial ©■ Note that if we set M = eT ^ gX , then by 
(j3J) we have 



E 



M 



\Z M (U, 



\2k 



(G(k + l)f 



N 



G(2k + 1) \e~t\ogX t 

which is the same answer we find in Theorem @J This is easily explained by the fact that in our 
model the eigenangles are multiplied by e 7 log X and so their mean density is M/2tt. Given that 
for random matrices the mean density is the only parameter in the asymptotics of local eigenvalue 
statistics, it is natural that the result should be the same as for unitary matrices of dimension M, 
since their eigenangles have precisely this mean density. 



Proof. Heine's identity j2J evaluates the expected value in (|21|) as a Toeplitz determinant 



E 



N 



' N 

n 

,n=l 



det 



(22) 



where 



1 

2^ 



<j>{9)e 



-m6 



d6 



is the nth Fourier coefficient of <p(9). The Toeplitz symbol (f)(9) is singular since it is zero when 
9 = 0. Thus, the asymptotic evaluation of this determinant requires knowledge of the Fisher- 
Hartwig Conjecture in a form proved by Basor pQ. 

We factor out the singularity in (f)(9) by writing 

(f)(9) = b(6) (2 -2 cos 9) k , 
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where 



6(0) =exp f-Hog(2- 2cos0) + 2k u(y) ( ^ Ci(|0 + 2vrj| logy log X)j dyj . (23) 

As we will see in the proof of Lemma below, the logarithmic singularities in the exponent on the 
right cancel. Thus 6(0) never equals zero. The asymptotic behavior of the Toeplitz determinant 
with these symbols has been determined by Basor pQ. She showed that if k > —1/2, then 

det [b-j]^^ ~ E exp £ log 6(0) de\ N k * (24) 



as N — > oo, where the constant E is given by 



\o g b{e)e- m6 &e) 1 6(o)- fcG2(A:+ " 



G(2A; + 1) 



To evaluate i£ we need to know 6(0) and the Fourier coefficients of log 6(0). These are given by 
the next two lemmas. 



Lemma 5. Let 6(0) be given by (J23J). Then 

— f log6(0)e" m9 d0 

271 " J -TV 



1 r,..^.-^,/, JO i/n = 0, 



w/iere 

/oo 
u(y) dy . 

Lemma 6. Let 6(0) 6e given by (|23|) and Zei u(x) Ztawe totaZ mass one wfft support in [e 1 ~ 1 ' x ,e]. 
Then 

6(0) = exp (2k (log log X + j)j(l + O k (j-L- 

Before proving the lemmas, we complete the proof of Theorem |1J Since u is a nonnegative 
function supported in [e l ~ l l x , e] of total mass one, we see that 



v(t) = 

and < v(t) < 1 if t € [e, e 1 " 1 ^]. Thus 



1 if e 1 - 1 ^ 
if t > e , 



n=l v v & // n<(l-l/A-)logX \(1- 



^ n 

l/X)logX<n<logX 



The first sum on the right equals log log X + 7 + O (1/logX) and the second is O (X x ) . Hence, 
we find that 

2 



1 ( n \ ( 1 

Using this and the value of 6(0) given by Lemma El we obtain 

E = exp (-*> glog * + 7) ) ^±ffi (, + 0l . (jJj)) 
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The proof of Theorem is completed by combining this, the case n = of Lemma 03 (1221) . and 
21. 



Proof of Lemma 03 We wish to evaluate 

^ f\ogb{e)e- m6 dO, 

where b{6) is given by (|23|) . After some straightforward algebra we see that this equals 

—k f w 2k f e ( f°° \ 

— J Iog(2-2cos0)cosn0d0 + — J u(y)lj Ci(01ogy log X) cos n6 dO J dy . (25) 

When n = the first integral vanishes by symmetry, and the second vanishes because 

POO 

/ Ci(0) d9 = . 
Jo 

This is a special case of the formula (see Gradshteyn and Ryzhik ^7 , p. 645) 



/>oo 

/ Ci(A0)cosn0 dO = < 
Jo 



-§i XA<n, 

-t ifA = n, (26) 
otherwise 



for A > 0, which we require below as well. Thus, both terms in (|25|) vanish and Lemma |S] holds in 
this case. 



When n is a positive integer, the first term in (|25[) equals 

for* k f n ( \ 

— / log (2 - 2 cos 6*) cos n0 d6> = / I log 4 + 21og(sin -) ) cos nO dO 

^ Jo n Jo V 2 / 

4A- r 77 / 2 

= / Iog(sin0)cos2n0d0 (27) 

n 

(see Gradshteyn and Ryzhik 7 , p. 584). The second term in (|23|l is, by (|26|). 

2k f e f f°° \ k f en/logX 

— / M (y) I / Ci(0 logy log X) cos n9 dO I dy = / u(y) dy 

k J i \Jo J nJi 



£ (v{e n ' l ^ x )-l) . (28) 



Inserting (|28j) and (|27jl into (|25j) . we find that for n > an integer, 

— \ogb{e)e- me d6 = - v (e n/losX ) . 
2vr J_ T n 

This completes the proof of Lemma El 

Proof of Lemma HO We calculate 6(0), where 



b{9) =exp ^-/clog(2-2cos0) + 2£;^' u(y) ( ^ Ci(|0 + 2vrj| logy log X)j dyj . (29) 



Using the expansion 

Ci(x) = 7 + log x + 0(x 2 ) 
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for x > 0, we find that the first term in the exponent and the j = term combined contribute 

- k log(2 - 2 cos 6) + 2k J u{y) Ci(|0| log y log X) dy 

= 2fc{-log(|0|) + O(# 2 ) + ^ u(y){\og{\e\\ogy\ogX)+ 1 + Ox(e 2 )) dy} 

= 2k | 7 + log log X + jf u(y) log logy dy + O x (# 2 

since has total mass one. Moreover, n(x) is supported in [e 1-1 /^, e], so we have 

1 



u{y) log log y dy < 
j A 



Therefore we find that 



Jim | -/clog(2-2cos0) + 2A; J u(y) Ci(|6>| logy log X) dy 



(30) 



2£(loglogX + 7 ) + O fe ( Y ) • 



Now consider the contribution of the terms with j ^ in (|29|) . An integration by parts shows 
that 



Ci(x) 



cost , sinx n . 1 . 
-— di = + — 

t x x A 



for a; positive and S> 1. Thus, since (1 — l/X)logX ^ logy log A" ^ logX, X > 2, and £ (— 7r,7r] 
we see that 



Ci (|0 + 2tt/| logy log X) 



1 



sin (|6> + 2ttj\ logy log X) ^ 



j=—oo 



i 



logy log X |0 + 27T7'| ' V (logX) 2 



In a standard way (via Abel partial summation), one can show that the series on the right is 
uniformly convergent for y € [e 1_1 / x , e], except possibly in the neighborhood of a finite number of 
points, and boundedly convergent over the whole interval. Moreover, the series may be bounded 
independently of £ (— 7r, 7r]. We may therefore multiply by the continuous function u(y) and 
integrate to find that 



\ 



Ci(\e + 27rj\ logy log X) 



j=-oo 



dy < 



1 



/ 



log X J 1 log y 



dy + O 



(logX) 2 / logX' 



uniformly for 9 £ (— it, it]. Combining this and (|30l) with (|29|) . wc obtain 

1 



6(0) = exp 2k (log log X + 7 ) + O k 



logX 



exp 2fc(loglogA + 7 ) l + O fe 



1 



logX 



This completes the proof of Lemma El 
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5. The Proof of Theorem 3 



First we prove Theorem when k = 1. In this case G 2 (k + 1)/G(2k + 1) = G 2 (2)/G(3) = 1, 
and by Lemma |U we may replace Px{\ + by Px(h + Thus, it suffices to show that for 
X«(logT) 2 -% 

As in the proof of Lemma El we write <S(X) = {n : p | n ==> p ^ X} and 

neS(X) 



where a_i(n) = /Lt(ra), the Mobius function, if n € S{yX); a_i(p) = for all p ^ X; and 
a_i(n) <C d 3 / 2 (n) <C d(n) for all n 6 <5(X). By (|16|) . if the e above is sufficiently small, we find 
that 

nG-S(A-) 



(The exponent 1/10 in place of 1/2 is accounted for by the slight difference between the conditions 
X <C (logT) 2_e and X <C (logT) 1 ^ 1 / 2+e \) Now for m and n coprime positive integers, we have 
the formula 



/ |C(i + ii)| 2 (-Y dt = -L= flog f— ^— ^ + 2 7 - + O (mnT 8 / 9 (log T) ( 
Jt \ n / \Jmn \ \lixmn J J \ 

(For example, see Corollary 24.5 of JJ.) Using this and the main term in (|31j) with 9 = 1/20, we 
find that 



i /-2T 



2 

E Q-i(n) 
n l/2+it 

neS(X) 



dt 



- ^ ~ ; K») log ; ^ +2 7 -l + -T-^logT) 6 , 

m,n£S(X) 

(32) 

where (m, n) denotes the greatest common divisor of m and n. The O-term contributes 

« r-v^iogr) 6 1 d(n) I « ^-^(logr) 8 . 
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Grouping together those m and n for which (m, n) = g, replacing m by gm and n by gn, and then 
using the inequality d(ab) ^ d(a)d(b), we find that 



E 

m,ndS(X) 



a_i(m) a_i(n) 



m n 



(m, n) I lo, 



(m, n) 2 
27rmn 



+ 2 7 - 1 



« E J E 

(m,n)=l 



d{gm)d{gn) log ran 



ran 



« E ^1 E 



c2(n) logra 



K n&S(X) 



n 



If we write f{a) = z2neS(X) d(n)n~ a = f\ P ^x (1 ~~ P~ a ) > then the sum over n is — / (1), which, 
by logarithmic differentiation, is 2/(1) J2 P ^x ^°SP/(p~ 1) ^ 7(l)(l°g^0 ^ (log^O 3 - We also have 
J2 g eS(X) d{g) 2 g~ l -C llp^x (l — 4 ^ (l°g^) 4 > an< i so the expression above is <C (logX) 10 . 

Thus far then, we have 

2 

a_i(n) 



1 r 2T 



n l/2+it 



df 



(33) 



logT ^ 



a_i(m) a_i(n) 



[m,n) 



m,n&S{X) 



m n 



+ 0((logX) 10 ) . 



Since Yl g \ n 0(d) = n > * ne remaining sum here is 



E 



a_i(m) a_i(n) 



m,n£iS(X) 



m n 



/ \ 



E 



0(5) 



\ff|n / <?e-S(X) \ neS(X) 

We wish to extend the sums on the right to all of S(X). For this we use several estimates. First, 



\ 



E 



a-x{gn) 



n < T V20„-l 



n 



(34) 



/ 



E 

ne-S(X) 



l"-i(s™)l 



n 



d(g) £ ^ = II f 1 - -) 2 «%)0og^) 



Second, 



E 



\a-i(gn)\ 



n>T i/2 0g -i 

neS(X) 



n 



< d{g) 



d{n) 



n>T l/20 9 -l 
n£S(X) 



n 



^d(g) 



/ j-l/20 " 



-1/4 



E 

n&S(X) 



d(n) 



n 



3/4 



« d(g)g^T-^ [J f 1 " 4i) _2 « %)5 1/4 T- 1 / 80 e 10 ^ 1/4 /l°^ 



say. From these it follows that the square of the sum over n in (|34|) is 



E 

,neS(X) 



a-i(gn) 



+ o(%)V /2 t- 1 / 200 



(35) 
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By arguments similar to those above we also find that 

<t>{g)d{g) 2 _ ^i/4oo , HaMg) 2 



E 

geS(X) 



9 



3/2 



«T 1 / 400 and Yl 2 «^ 1/10 °- 



9>T l/20 

g£S(X) 



g- 



Using these and (J22J), we find that the right-hand side of equals 

/ \ 



E - E 

g&S{X) 9> t 1 /20 

V geS(X)J 



4>{g) ( a-i(gn) j + () 



( \ 

4>{g)d{g? 



\neS(X) 



Ef IE 



2^-1/200 ^ 



V 



96<S(X) 



,3/2 



/ 



a-i(pn) 



9 £5(A) y VneSpC) 



Combining this with (j33|) . we now have 



22" 



T 



|C(| + i*)| 



E Q-i(n) 
n l/2+ii 



di 



logT 



E 



a_i(gn) 



gdS{X) * \naS{X) 



+ 0((logX) 10 ) . 
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(36) 



Since a_i and are multiplicative functions, we may expand the entire sum into the Euler product 



LT EEE 

p^X \ r j k 



P 



2r+j+k 



Recall that a_i(n) = fJ-(n), the Mobius function, if n G S(yX); a_i(p) = fj,(p) for all p ^ X; and 
a_i(n) <C d 3 / 2 (w) <C d(n) for all n G Thus, the product equals 

nK) n .K+°(?))-nH) n 



v / X<P<X 



l 



eTlogX 



l + O 



1 



logX 



Since log X <C log log T, it now follows from ()36|) that 

2 



1 



2T 



E 

neS(X) 



a_i(n) 



1/2+it 



eT log X 



logX 



(37) 
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Rewriting (jSH (with 9 = 1/20) as P* x {\ + it)" 1 = E + 0{T- e / 200 ), we see that 

-e/200/° 



1 

T 



,-2T i p2T , 

ic(i+it) J pi(i+itr i i 2 d*=- ic(*+it)i 2 |E+°( r 



1 f 2T | 1 2 /I /•» , , \ 

+ °(^7M/ T 2T K(i + ^)l 2 dt). 

The final term is O(T -<E / 20 °) since the second moment of the zeta function is 0(T log T). Also, by 
the Cauchy-Schwarz inequality and (|37j). the second term is 

1 / f 2T 1 2 r 2T \ 1 ^ 2 

«j^m>{j T |c(i+*)E| d */ T lc(l+^)| 2 dtj 

T 2 iog 2 r/io g x) 1/2 «T- £ / 400 . 



1 



^l+e/200 

From these estimates and Q37[l. we may now conclude that 

i /; t i«i + + * = ^ (i + o )) 

for X = 0(logT) 2_e . This completes the proof of Theorem 01 in the case k = 1. 

We now prove Theorem for k = 2. By Lemma 121 we may again replace P x ( 3 + it) by P£ (3 + it) . 
Furthermore, G 2 (3)/G(5) = 1/12, so it suffices to show that 



for X <C (logT) 2_e . By ()16|) (see (|31j) also and the remark following it), we have 



(38) 



E 

neS(X) 



n 



+ O ( r - e e/io 



say, where a_2(f>) = —2 for all p ^ X, a_2(p 2 ) = 1 if p ^ VA, a-2(p 2 ) = 2 if \fX < p ^ X, and 
a_2(p : ') = otherwise. In particular, we note that |a_2(w)| ^ d(n). 

In carrying out the proof of splitting for this case, we will gloss over some of the less important 
steps as these are handled analogously to those for the k = 1 case. In particular, by an argument 
similar to the one at the end of the proof of the case k = 1, one can show that 

i-2T 



1 r zl 



2 D * (l 
2 



it 



l + O 



1 



logX 



,-2l 2 



At 



2T 



(39) 



E 



<x_ 2 (n) 



neS(X) 



11 



1/2+it 



df . 



where Y = T d and > 0. Eventually we will take 9 very small. 

To estimate the right-hand side we use an analogue of (|32|) due to Jose Gaggero [12 . Let 
A(s) = J2 n <Y a n^ s 5 where the a n are complex coefficients and Y = T e with 9 < 1/150. Gaggero's 
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formula is 
/ 1 \ 1 f' 2T 

4 

=E 



k=i 



c k (m,n)a m a n , k (YT(m,n)\ k fT(m,n) 

m,n<Y x 7 x 



_ £ _ M^LL + (i)) V- V i 

^ mn ^ d V 4eT y ' j ^ v ^ u 

m,n^Y 0<d<Y/4 K ' v<Vi u<Ui 

(n d u,m d v)=l 

Ea m a^ \ - (m,d)(n,d) / mn . \ ^ 1 x - 1 

^ d l lo s W + 0(1) J ^ - 2. u- 

m,n<y 0<d<mn/4Y ixrV/ u<f/{ 

(n d u,m d t>)=l 

Here 

Ui = CYT/dn d , Vi = CYT/dm d: (40) 
C/{ = CmnT/Ydn d , V[ = CmnT/Ydm d , 

C = 2/tt, B is an arbitrary positive number, and for integers n and d we write n d = n/(n, d). Also, 
Ci{m,n) = (1/Air 2 )5(m n )5(n m ), where 



IjA (1+1/P)J 



and Cj(m,n) <C |c4(m, n)|(loglog3mra) 4 ■? for j = 1,2,3. 

To estimate the right-hand side of (j39j) . we take a n = a_2(n) and y = T ei in this and obtain 



( 1 + (six))fX IC(U^(§ + . t )-'f * 

= (^+£X«))vt e 



a_2("i)o-2(^)<5(?ri/ (ret, n))5(n/ (m, re)) . . 

{ T7Z . 77/ J 



m,ne<S(A") 

_ a_ 2 (m)a_ 2 (n) ^ (m,d)(n,d) / F + \ 1 1 

m.n^Y 0<<2<Y/4 v 7 u<V a u<E/i 

m,ne5(A) (n d u,m d v)=l 

\- a_ 2 (ra)a_ 2 (n) (m,d)(n,d) / ,mn , ,\ 1 \- 1 

m,n^y 0<d<mn/AY v<V{ u<U[ 

m,n£S(X) (n d u,m d v)=l 

Tx-T 2 - r 3 , 



say. 
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Let us denote the sum in Ty by Si- Grouping together those terms for which (m,n) 
then replacing m by mg and n by ng, we obtain 



g and 



/ 



h= E ] E 

gGS(X) n&S(X) 



a- 2 (gn)5 (n) 



11 



\ 



E 



m^Y/g 
(m,n)=l 
\me5(X) 



(42) 



Let P = Ilp^xP- Since a_2 is supported on cube-free integers, the g's we are summing over 
may be restricted to numbers of the form 

9 = 91 dli where gi \ P, g 2 \ (P/fli) • 

Note that this representation is unique and that (51,(72) = 1- The summation over g in (|42j) may 
therefore be replaced by the double sum 

E E • 

911 92|(P/Sl) 

In the sum over n we group terms together according to their greatest common divisor with 
9 = 9x92' Observe that we may assume that (n, g 2 ) = 1, for otherwise a cube divides gig\n and 
a -2(9 n ) vanishes. If we then write (n, g{) = r and n = rN, we may replace the sum over n in (|42|) 

by 

E E ■ 

r\9i N^(Y/ rgi g%) 
NES(X) 
(N,(gi/r) 9 2)=X 

Ignoring the restriction (m, n) = 1 for the moment, we may similarly write the sum over m in 
flUD as 

E E • 

s\gi M^(Y/sg ig %) 
MeS(X) 
(M,(gi/s)g 2 )=X 

Instead of (m, n) = 1 we now have (sM, rN) = 1 or, equivalently, (M, N) = (r, s) = (N, s) = 
(M,r) = 1. We may impose the condition (r,s) = 1 by replacing s \ gi in Q by s | (gi/r) since 51 
is square-free. Furthermore, since (N,gi/r) = 1 and s \ (gi/r), we automatically have (N,s) = 1. 
Thus, the coprimality conditions on M are (M, (ffi/s)^) = (M,N) = (M,r) = 1. The first 
condition implies the third because r | [gi/s). Thus, we need only require that (M, N(gi/s)g<2) = 1. 
The sum over m may therefore be written 

E E 

s\ 9l /r M^Y/sgrgl) 
M&S{X) 
(M,N( gi /s)g 2 )=X 
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We now have 
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S 



9l fp ^07^ r|si s|(9l/r) 



«_ 2 (H^iV( 5l /r))(5(r7V) 



S2|(P/3l) 



NeS(x) 

(N,(gi/r)g 2 )=l 



E 



a_ 2 ( S 2 52 2 M( 5l / S ))5( S M) 



M<(y/s 9l9 2) 
MeSpO 
(M,Ar( Sl / s ) S2 )=i 



M 



Note that if N and r have a common factor, then a_ 2 (r 2 g%N (g\ / r)) = 0, and similarly for M and 
s. We may therefore replace the coprimality conditions in the sums over N and M by (N,g\g 2 ) = 1 
and (M, Ngig 2 ) = 1, respectively. The new conditions then imply that a_ 2 (r 2 g 2 iV(gi/r)) = 
a_ 2 (r 2 )a_ 2 (g2)a-2(^)a-2(si/ r '), S(rN) = S(r)S(N), and similarly for a- 2 {b 2 g\M{g\j ' s)) and 5(sM). 
Hence 



«- 2 (<7i) 2 v- a-2(52 2 ) 2 v-«-2(r 2 )5(r) V- a- 2 (s 2 )<*00 



31^ 



51 



E 



(iV, S i S2 )=l 



E 



<?2<(Y/ 9 i)2 
92\(P/gi) 

a^ 2 (N)S(N) 

N 



9l 



E «-2 
a. 



E 



-2(r)r f-; a_ 2 (s)s 

r\gi s\(gi/r) 



E 



a_ 2 (M)5(M) 



M<:(Y/s 919 2) 

(M,Ng l92 )=l 



M 



We next extend each of the sums here to all of S(X). The error terms this introduces are handled 
as they were in the case k = 1, and they contribute at most "little o" of the main term. Observing 
also that M and N may be restricted to cube-free integers, we obtain 

Sl = (1+o(1))E ^0! E ^ e ^im E m 



gi\P 



91 



E 

N\{P/gig2? 



92|(P/fll) 52 

a- 2 (N)6(N) 



a-2(r)r f-; a_ 2 (s)s 



TV 



E 

M|(P/JV 9 i 92 ) 2 



a_ 2 (M)(5(M) 
M 



We now define the following multiplicative functions: 

^ a-2(d)5(d) = Y[ (l + a -2(p)*(p) 

a_ 2 (d)5(d) 



B(n 
C(n 
D(n 
E{n 



p"-\\n 



E 



<L4(ci 2 ) 



n i+ 



d\n p"||n 

a_ 2 ((i 2 )(5(d) 



P 

a- 2 (p)6(p) 



■■■ + 



a- 2 (p a )6(p a ) 

p a 



pA{p 2 



+ ■■■ + 



a- 2 (p a )5(p a ] 



E 

d\(n,P) 



a^ 2 (d)d 
a_ 2 {d 2 )5{d) 



n 

p|(n,P) 



1 + 



a_ 2 



a- 2 (d)C(d)d H I 

d|(n,P) 2V 7 V 7 p|(n,P) 



a- 2 (p)p J ' 

Q„ 2 (p 2 )(5(p) \ 



E 

<i|n 



at- 2 (d 



>2\2 



A{d 2 )B{d?)d? 



P \l 



n (i+ 



a- 2 (p)C(p)pJ ' 
a-2(p 2 ) 2 



^(p 2 )5(p 2 )p 2 
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and 

v a- 2 (d) 2 C(d)D(d) _ rr / a_ 2 (p) 2 C(p)D(p) 
{ ' ^ A(d 2 )B(d 2 )E(d)d Hi A{p 2 )B{p 2 )E{p)d 

d\n p\l 

Using these definitions and working from the inside out in (|43|) . we find first that the sum over 
M is A((P / N gig 2 ) 2 ) = A(P 2 ) / A{N 2 )A(g\)A{gl). The contribution of the sums over M and N 
together is then (A(P 2 )/A(g 2 )A(g%)) (B{P 2 )/B(gj)B{g%)). Thus, so far we have 

a- 2 {gi) 2 a^ 2 (gj) 2 



\ - a- 2 (r 2 )^(r) a_ 2 (s 2 )<5(s) 



a_ 2 (r)r ^ a_ 2 (s)s 

The sums over r and s contribute C(gi)D(gi), and the sum over g 2 is then E{P) / E{g\). Thus, we 
see that 



9i\P 

?2 



q_ 2 (ff 1 ) 2 C(g 1 )ZJ(g 1 ) 
gi A(g 2 )B(g 2 )E( gi ) 



= (1 + o(l)) A(P Z )B(P Z )E(P)F(P) . 
Using the expression for F(P) as a product, we see that this is the same as 

* = (i + o(i)) n ( mpW)e( P) + o^mMm ) . m 

By the definitions of C and D we see that 



C{p)D(p) = 11 + + 



_ -, _ a- 2 (p 2 )^(p) 

since a_ 2 (p) = —2 for p dividing P. Similarly, 

«- 2 (p 2 ) 2 



A(p z ) B(p z ) E(p) = A(p z ) B(p') + 



It is clear that A(p 2 ) = A(p s ) = ■ ■ ■ . Therefore 



B( p 2 ) = 1 + a ~^P) 5 (p) + a- 2 (p 2 )5(p 2 ) 



Aip 2 ) 
and 



pA(p 2 ) p 2 A{p 2 

1 

2 



A{p 2 ) B(p 2 ) E{p) = 2A(p 2 ) - 1 + a ~ 2i f )2 



(45) 
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We use this, (|45[). and a_ 2 (p) = —2, and obtain 



^- 1 + i^£ + l.^W) 

_ 2%) + g„ 2 (p 2 )£(p 2 ) \ _ i + g_ 2 (p 2 ) 2 + 4 4a„ 2 (p 2 ) 2 5(p) 



p p 2 J p 2 p p 2 

^ i 4 - U(p) a- 2 (p 2 ) («- 2 (p 2 ) ~ 4a_ 2 (p 2 )J(p) + 2,5(p 2 )) 
p p 2 

Recall that 5(p r ) = l + r j}~^j , so that 5(p) = 2/(1 + 1/p) and <5(p 2 ) = 25{p) - 1. Also recall that 
a-2(p 2 ) = 1 if P ^ Thus, for p ^ Va the last line is 

= 1 + 4- 45(p) | 1 - 4g(g) + 2 (2£(p) - l) _ i | 4 - 4J(p) 1 

p p 2 p p2 

= 1 + i 8 1 _ (1 - 1/p) 3 _ (1 - 1/p) 4 



p p+1 p 2 1 + 1/p 1 — 1/p 2 
On the other hand, if \fX < p ^ X, then a_ 2 (p 2 ) = 2, and the last line is 

_ (1 - 
1 - 1/p 2 

Combining these results in (|44jl. we find that 



4 



+ O (VP 2 ) • 



*-a + *))n(^) n 
=( i +°( i »n( I ^f ! ) n (i+o (VP 2 )) 



Since 



(l+otiDn^-i/pfnii-v/r 1 

P^X p 

(l + (l))^(e^ogA)- 4 . 
b 



Ti = s 1 (JL + o( ei ))io S J :/' 



and ei > may be taken as small as we like, we now see that 



1 , A / logT x 1 



To treat the second term on the right-hand side of 1)41(1 . 7^, we require two lemmas. 
Lemma 7. Suppose that a and b are positive integers with (a, b) = 1. Then for b ^ x, we have 

i _ m 



V - = ^Vlogx + 0(loglog26). 
^— ' n b 

(an,b)=l 
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Proof. Since (a, 6) = 1, the condition (an, 6) = 1 is equivalent to (n, b) = 1. Thus, the sum is 

n^a: n^x d|n d|b m^x/d d\b 

(n,b)=l d\b 



Now EJ d | b u((i)/d = <j)(b)/b and 



u(d)logd _ 0(6) ^ logp ^ 0(6) 

Furthermore, 



V 7 , = ^ V —^7 < ^ log log 26 « log log 26. 
a b ^— ' p — 1 6 

d|6 p|6 



£^nK)^<<~- 



Thus we find that 



d\b p\b 



V i = ^logx + 0(loglog26). 
n 6 

(an,ft)=l 



-1 



□ 



Lemma 8. Zei k(ji) = Y[ p \ n ( 1 + |j and* ~ei [/, V 6e either U\, V\ or U[, V{ as defined in (|40[). 
If m,n,d <^.Y t 11 / 150 ^ a nd (n^, m^) = 1, i/ien 

S - = 4 K ( m ^) K ( n rf) 1 °g ?71o g T/ + ( lo g Tlo g lo g T )- 

(n d M,m ti t))=l 

Proof. The conditions (n^it, m^) = 1 and {rid, rrid) = 1 are equivalent to (v, rid) = 1 an d (it, rridv) = 
1. Hence, by Lemma [7] the double sum equals 

£ - £ ± = iogtf £ if^l + ooogiogK^ 

v „ u *—? r v V rndV 

v<V u<U v<V v ' 

(u,n<j)=l (u,m d v)=l (v,n d )=l 

= \o S U £ -(^^)+0(lo g Tloglo g T). (47) 
~(, v \ m d v J 

(v,n d )=l 

Denoting the sum on the right by ">""", we have 

v<V r\m d v r<m d V v<V 

(v,n d )=l ( v ,n d )=l 

r\m d v 

Now set {rrid, r) = g and write r = gR. Then {rrid/g, R) = 1 and we find that 

^ q 2-~< R ^ v' 

g\m d a R<m d V/g v<V 

(m d /g,R)=l W"d)=l 
R\v 
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If we set v = Rw, then w <V/R, and (Rw,n d ) = 1 is the same as the two conditions (R, n d ) = 1 
and (w,rid) = 1. Thus, using Lemma and the observation that the inner sum vanishes unless 
R < V, we obtain 



g\m d * R<m d V/g w<V/R 



(m d /g,R)=l 
(n d ,R)=l 



E 1 

<V/R 
(w,n d )=l 



g|m d y R<U \ a / 

(m d /g,R)=l 
{n d ,R)=l 

We may assume (R,g) = 1, for otherwise fi(gR) = 0. The coprimality conditions on the sum may 
then be written (m^n^, i?) = 1, and we find that 



E = E^ E ^(^4+o mnd) 

g\m d y R<V V a 

(R,m d n d )=l 



4>( n d) 

n d 



(R,m d ?i d )=l 



+ f.oglo g 2n i y:^E i! l#) 

\ 9 |m d y R<V J 



Since E 9 |m d KsOI/fi' = H p \ md ( 1 + 1 /p) < foglog2m d , the error term is <C (loglog2m d loglog2n d ). 
The main term is 



rid 



g\m d 



E 



, ij=i 

\(R,m d n d )=l 



R 2 



+ o(v- 1 ) 



J 



cor 1 n (^V^^E^ i^ E i^)i 

P\m d n d g\m d \ g\m a 



V ^ g 

g\m d 



7T 



n H 



P\m d n d 



[m d ) <p{n d ) f'logV log log 2m d 

log 1/ + C 



By hypothesis, (m d ,n d ) = 1. Furthermore, J] p |i (l - Vp 2 ) 1 = U P \i i 1 + Vp) 1 = K (0- 

Thus, combining our estimates, we obtain 



y~) = -^K(m d )K(n d ) log V + 0(log log 2m d log log 2n d ) . 

Since m, n <C T 1 / 150 , we obtain from this and (|47|) that 

1 16 
V- V - = -zii{m d )K{n d ) log logy + 0(logT log log T) . 

Z J 71 Z S 71 tt^ 



f ' U IT' 

v<V u<U 

(n d u,m d v)=l 



□ 
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Returning to 7^ in (|41|) and using Lemma |HJ we have 

6 y «- 2 (m)a- 2 (n) y (m,d)(n,d) / Y 
7r 2 ^ mn ^ d V 4cr 

m,n^Y 0<d<y/4 v 

m,ne5(X) deS(X) 

x U(m d ) K (re d ) log [/1 log Fi + 0(log T log log T) 

where XJ\ = CYTjdn^Vx = CYT/dm d , and Y = Tj\ Interchanging the order of summation, we 
find that 



* = £ E K log( s> +o(1 E 



a_2(m)a_2(n)(m, d)(n, d) 



dES(X) m,n&S(X) 

x f «(m d ) K (n d ) log C/i log Vi + 0(log T log log T) 

Since k(ti) <C log log 3n, the expression in the last parentheses is 

= K {m d ) K {n d ) (logtfi logVl+0(logrioglog 3 T)) = (1 + 0(ei))K(m d ) /c(ra d ) log 2 T. 

Thus, 

2 



r 2 = 4(i + o (ei ))io g 2 r £ "° g(r/4rf i + 0(1>) 



7T 

o<d<y/4 



E 



ck-2(w) (n, d)n(n/ (re, d)) 



n 

. n<:y . 
\neS(X) J 



(48) 



Denote the inner sum by S(d). As on previous occasions, extending the sum to all of S(X), we 
introduce an error term that is o(l) times the main term. Thus, grouping together terms in S(d) 
for which (re, d) = e, say, we obtain 

S(d) = (l + o(l))5> £ a - 2(n) re K(n/e) =(l + o(l))^ £ a - 2(e ^ )K(iV) . 

e|d neS(X) " e|d NeS(X) 

(n,d)=e {N,d/e)=l 

Since a_2 is supported only on cube-free numbers in S(X), we may assume that e \ P 2 . Therefore, 
e I (d, P 2 ) = D, say. Now D may be written uniquely as D = D\D\, where D\ \ P and D2 \ {P/D%), 
so that, in particular, (D\,D2) = 1- Furthermore, we may write any divisor e of D as e = e^ejj, 
where ei | D\,e 2 \ D 2 , and es \ (D 2 /e 2 ). Note that this means the are pairwise coprime. The 
condition (N,d/e) = 1 is now (iV, (DiD^/e^e 2 )) = 1. Also, a_2(eiV) = a^^i^eliV), so we 
may assume that (N, e$) = 1 and, therefore, that (N, (D\D 2 / e\e 2 )) = 1- Observe, moreover, that 
e 2 I £> 2 implies e 2 | (^l/e 2 ). Thus, (iV, (Di-Df/e^)) = 1 is the same as (N,(D 1 D 2 /e 1 )) = 1. It 
follows that N and ei can have a common factor, but not N and 62 or 63. We may therefore write 
a-2(eie 2 e§iV") = a_ 2 (eiiV)a-_2(e2)a-2(el) and 

5(d) = S(D) = (1 + 0(1)) £ £ a_2(e 2 ) £ a_ 2 (e|) £ «-*MO"W 

ei|D ie2 |D 2 e 3 |(D 2 /e 2 ) NeS(X) 

{N,{D 1 D 2 /e 1 ))=l 

= (l + o(l)) 2^ ]v ^ «-2( e 3) 2^ a -a(ea)- 

ei|Di 7V&S(X) ' e 3 |D 2 e 2 |(D 2 /e 3 ) 

(JV,(£) 1 D 2 /ei))=l 
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The innermost sum is 

e «-2te)= n (i+a- 2 (p))= n (1-2) 

e 2 |p2/e 3 ) P|(0 2 /e 3 ) p|(£> 2 /e 3 ) 

= fJ,(D 2 /e 3 ) = V(D 2 )n(e 3 ). 

We also have 

E Me 3 )a-2(e|) = [I " ^(p 2 )) . 

e3|£>2 p|£>2 

At this point it is convenient to define numbers 

Pi= U P and P 2 = 11 p. 

Notice that P = P\P 2 . Since ot- 2 (p 2 ) = 1 if p \ Pi and a^ 2 (p 2 ) = 2 if p | P 2 , the sum over e 3 
equals unless D 2 \ P 2 , in which case it equals n(D 2 ). Thus, if D 2 and Pi have a common factor, 
5(Di£)|) = 0, whereas if D 2 | P2, then 

flW = S ( Dl Di) = (i+o(i)) yj £ a -^y K(iV) ■ 

ei|£»i NeS(X) 

(AT,(DiD 2 /ei))=l 

From this point on we shall therefore assume that D 2 \ P 2 . 
Now set (N, ei) = r and write ./V = rM. Then we have 

s( z> lD 5) = (i + °(D)££ £ °- 2(ei ^ ) " (A ' ) 

ei|Dir|ei ATeS(X) 
(JV,ei)=r 
(AT,(D 1 D 2 /ei))=l 

= (i+ (i))EE- E m • 

ei|Di r|ei MeSpO 
(M,ei/r)=l 
(rAf,(Di£> 2 /ei))=l 

We may assume that (M, r) = 1 and (r,e±/r) = 1, since otherwise a_2(r 2 M(ei/r)) = 0. Actually, 
(r, ei/r) = 1 is automatically satisfied because r | e\ and e\ is square- free. It follows that n(rM) = 
n(r)n(M) and, since we also have (M, e\/r) = 1, that a_2(r 2 M(ei/r)) = a_2( r2 ) a -2(-^0 a -2( e i/ r )- 
The coprimality conditions in the sum are now seen to be equivalent to the conditions (M, r) = 
(r, ei/r) = (M, e\jr) = (M, (DiD 2 /e\)) = (r, (DiD 2 /e\)) = 1. As we have already pointed out, 
the second of these is automatic. Similarly, so is the last. The remaining conditions are equivalent 
to (M, DiD 2 ) = 1, so we find that 

Q(n rfl\ (1 , ji^r. )r a -2( r2 ) K ( r ) V" q- 2 (M)K(M) 
S(£>ii3 2 ) = (1 + o(l)) ^ a- 2 (ei) ^ (r) E m ' 

ei|Di r|ei V ; MeS(X) 

(M,DiD 2 )=l 



The sum over M equals 



II ( 1 + a - 2{p) p K{p) + Q - 2( ^ K(p2) ) = G(P/ Dl D 2 ) , (49) 



p|(P/DiD 2 ) 
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say. Hence, 

S{D x Dl) = (1 + o(l)) G(P/D 1 D 2 ) £ a_ 2 ( ei ) £ "-^"fr) 

zi\Dx r[ei 2V ; 

The double sum equals 

E a ' 2(r2) A (r> E ^(/ 1 ) = MP.)E Mr) "' 2 M' c(r) 



_ MA)n (l + f^g=W) (50) 
= n(Di)H(Di) , 

say. Thus, 

S(<i ) = S ( DlD 5) = (l + (l,)G(P)^™, (51) 

provided Z?2 | P%\ otherwise S(d) = 0. 

We use this in (gSJ). Recall that for each d < Y/4 we had set (d, P 2 ) = D\D\ with D\ \ P and 
D 2 | (P/Di). Recall also that P = P X P 2 and K = T ei . We therefore have that 

r 2 = 4(i + 0(ei ))iog 2 T E ( '° g(F/4rf J + 0(1)> 5(^ 

0<d<y/4 

= ^( 1 + 0(£l))log2rE _L E « E WW^D+od)) 

£> 2 |P 2 2 Dx\{P/D 2 ) 1 0<5<y/4DiZ)2 

(5,(PiP 2 ) 2 /£ , i^)=l 

The coprimality condition in the last sum is equivalent to (5, PiP 2 / D 2 ) = 1. Thus, using (foil) , we 
find that 

^'(HO^fV^^ E 

P> 2 |P 2 2 V 2; D 1 \(P/D 2 ) V 17 

(log(y/(4 J D lJ D 2 5)) + 0(l)) 



x 



E 

0<5<y/4DiD| 
(5,PiP 2 /D 2 )=l 



By Lemma [7] the sum over 5 is 



1 0(P) L> 2 9 

« logy e j«V^ logy ' 



{8,P 1 P 2 /D 2 )=l 

Thus 



T 2 « G (P) 2 ^i og2 rio g 2 y E E 

P> 2 |P 2 ^ 7 V ' D t \{P/D 2 ) V 



If we denote the innermost sum by I(P/D 2 ), then 

pG(p) 
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and we find that 

Now, by the definitions of G, H, and / in (JUJ), (j5U|). and ([52]). we have = 1 - | + O(^), 
ff(p) = 1 + 0(1), and I(p) = 1 + I ((1 + 0(|))/(1 - | + O(^))) 2 = 1 + J + 0(£). From these 
estimates it is clear that the product over p dividing P 2 here is Uvx< P ^X i 1 + (Vp 2 )) < 1- Thus 

T 2 « e ?G(P) 2 /(P)^log 4 T 



« 6?log 4 Tn ((l " - + 0(l/p 2 )) (1 + 0(1//))) 
p\p P '' 1 

« £ ?io g *rn(i-£) 4 « ei 



log T N 4 



logX 



The treatment of T% is almost identical and leads to the same bound. Thus, combining our 
estimates for T\ (see (|4fi|)). I2, and T3 with (|41|). and noting that we may take ei > as small as 
we like, we obtain IfMHJI , This completes the proof of the case fc = 2 of Theorem |31 and thus, also 
the proof of the theorem. 



Appendix A. Graphs 

To illustrate Theorem^ in Figures ^-[3] we have plotted |Zx(|+rt)| and |-Px"(f fo r t near the 
10 12 th zero for two values of X, and have compared their product with the Riemann zeta function. 
The values of X used are X = 26.31 ~ log 71912 and X = 1000. Though the functions Px and Zx 
depend upon X, when multiplied together the X dependence mostly cancels out, and we have an 
accurate pointwise approximation to the zeta function. The actual functions plotted are 



\Px (| + i(x + t )) I = exp 




A(n) cos((x + to) log n ) 
log riy/n 



and 



N+10Q 

\Z x {\+i(x + t Q ))\= Yl exp(Ci(|x + t -7n|logA)) , 

n=N+l 

where to = 7io 12 +40- Fhe values of the zeros of the zeta function came from Andrew Odlyzko's 
tables |20| . The functions were plotted for x between and 5, a range covering the zeros between 
7io 12 +40 an d 7io 12 +60- Note that the function Zx we have plotted is an unsmoothed, truncated 
form of the function Zx that appears in Theorem ^ 
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1 2 3 4 5 



Figure 1. Graph of \((± + i(x + t ))\ (solid) and \P x (^+i{x+t ))Z x (^+i{x+t ))\, 
with to = 7io 12 +40i w ith X = logio (dots) and X = 1000 (dash-dots). 



1 2 3 4 5 



Figure 2. Graph of \Px{\ + i{x + io))|, with t = 7ioi2+4o> wit h X = logt (dots) 
and X = 1000 (dash-dots). 

6r 



fl l\ 

.Ann /U\a *M M;U 

12 3 



Figure 3. Graph of \Z X {\ + i(x + t ))\, with t = 7ioi2+4o> wit h X = logi (dots) 
and X = 1000 (dash-dots). 
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